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Anomalous transport and reaction dynamics are considered by providing the theoretical grounds 
for the possible experimental realization of actin polymerization in comb-like geometry. Two limiting 
regimes are recovered, depending on the concentration of reagents (magnesium and actin). These are 
both the failure of the reaction front propagation and a finite speed corresponding to the Fisher-KPP 
long time asymptotic regime. 
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I. INTRODUCTION 

Microfluidics is an indispensable tool of modern bio-physical research. It allows to perform complex single-cell 
experiments with an immense throughput and high level of control. A flexible design allows for custom geometries 
and control of flows and chemical reactions. Recently, to probe the dynamics of actin polymerization, as well as to 
use the geometry of microfluidic device having the main supply channel with numerous identical side channels or 
chambers of different shapes The main channel serves to deliver and fill the side chambers with reagents where the 
corresponding reacting can be and observed. The flow in the main channel and diffusion in the side-channels are 
dominating means of transport in such devices. Remarkably, the process of diffusion in this particular geometry was 
extensively studied in the context of anomalous diffusion. It is known as a comb model and it was demonstrated 
that the transport of particles along the main channel (called backbone in the model) can become subdiffusive when 
the particles get trapped by diffusing into the side channels. Until recently it was mostly an abstract model, which 
was, however, extremely useful in understanding the principles of anomalous subdiffusive transport (many references 
here). The experimental setup on actin polymerization [li, 2) is the direct implementation of the comb model 
where the effects of complex diffusion should have a substantial effect on the observed phenomena. Interestingly, the 
comb structure not only leads to an anomaly in transport but also to a very remarkable effects on the propagation of 
chemical reactions 

The goal of this paper is to combine the consideration of anomalous transport and reaction dynamics to provide 
the theoretical grounds for the corresponding experimental efforts. Our analytical results on reaction propagation can 
help to guide the design of microfluidic devices but also can lead to real experimental tests of anomalous diffusion and 
reaction dynamics. For the reaction of polymerization, depending on the concentration of reagents we can recover 
such remarkable phenomena as the failure of reaction front propagation |7, || or a finite speed or a Fisher-KPP long 
time asymptotic regime [§, Rj. 


II. MAPPING OF THE LAPLACE OPERATOR ON A COMB EQUATION 

Mapping of the Laplace operators in the 3 D complicated continuous-discrete geometry on a continuous comb 
model equation, relates to averaging over xyz- space. In reality, we have two Laplace operators acting on distribution 
functions in a bulk P b (x, y, z) and in fingers (where reactions take place) Pf(x, y, z). 

( i ) In the bulk of infinite length along the x coordinate and yz surface cross-section a x a we have for the Laplace 
operator 


DAP b (x, y , z) = D{&1 + d 2 y + d 2 )P b ( x, y , z) 
with the diffusivity of magnesium D and the boundary conditions 

d x P \ x =o = d x P\x=a = d y P\y=o = dyP\ y = a = 0, 


the following algorithm of mapping can be suggested. 
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Integration over z leads to the disappearance of the z component due to the boundary condition. Integration over 
y in the bulk yields zero except those y regions where the bulk is connected with the fingers. Plunging the fingers 
inside the bulk, one obtains that this region is at y = 0. Therefore, we have 

DAP b (x, y, z) => aDd 2 x P(x , y = 0). 

(ii) The fingers have height ft. (along the y) and the xz cross-section is of size ft x ft. Therefore, we have the following 
Laplace operator with diffusivity d inside the fingers 

dAP f (x,y,z) = d{d 2 + + d 2 )P f (x,y,z). 

Integration/averaging over x and z leads to zero for everywhere, except for d 2 in periodic regions 

dyP( x , y)^2S(x- nL) =>• ^-d 2 P{x, y ), 

n 

where L is a distance between the fingers. Combining both equations, one obtains a kinetic equation for the mag¬ 
nesium transport along the comb structure with parameters related to experimental geometry a, ft, L, ft and diffusion 
coefficients D and d 


d t P = DS(y/h)d 2 P + pdd 2 P. (1) 

Here p = bh/La is a dimensionless density of the fingers. In what follows we will use the notation for the coefficient 

d p = d ■ p. 

However, the scaling of the S function by ft can be inconvenient in the limit h —> oo. Therefore, in what necessary, 
we scale the S function by 

h b = Vhb , 6(y/h) 6(y/hb). (2) 

In this case, we have h —> oo and ft —> 0, and in these limits hb = const. This scaling supports the density p. Say 
hb = aL and in this subdiffusive limit p = 1. 

The boundary conditions are 

P(x = ±oo, y, t ) = d x P{x = ±oo, y,t) = 0 and d y P(x, y = ±h, t) = 0. 

The initial condition is 

P{x,y,t = 0) = P 0 (x)S(y). 

Amending comb equation m with reaction term inside fingers p ■ C(P), we have 

dtP = DS(y/h b )d 2 P ± d p d 2 y P - P C(P) , (3) 

Probably, another mapping on the comb model can be made in the framework of a Dykhne model dMl. In 

7 2 L 

this case, dimensionless density is p = ^prp- In some extend, this density sounds more physically. But this leads to 
additional argumentation of limits, for example ft. —> oo and b —> 0. Namely hb 2 —> a 2 L and h b —> h bb = [ftft 2 ]^. 

III. FEW POINTS ON REACTION FOR DEFINITION OF THE REACTION TERM 

For the reaction-polymerization, we use the following stoichiometry expression 

A ± B => C , 

where 

[A] is concentration of actin 

[B\ is concentration of magnesium 

[C\ is concentration of polymer 
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In general case, we consider the second order reaction, and also take into account that [B\ = P(x,y,t). Therefore, 
we have 


dP 

dt 


-kP ■ [A] 


with the the reaction rate k and initial condition [B(t = 0)] = P(t = 0) = Po, [A(t = 0)] = [A]o = N , [C]o = 0. 
Let us express [A] by P. We have 


P=[B\ = [B\ o - [C] = P 0 - [C], 
[A] = [A] 0 — [ C] = [A] 0 — Pq + P ■ 


Therefore the reaction equation reads 


dP 

dt 


= —k(N — P 0 + P)P. 


(4) 


If during the reaction time, N — Pq P and N Po, the second order reaction equation ((dj) becomes of the first 
order 


dP 

dt 


= -kNP. 


Substituting Eq. © in Eq. © and denoting C = p ■ k ■ N yields 

d t P = Dh b S(y)d 2 x P + d p d 2 y P - CP. 


(5) 


( 6 ) 


IV. SOLUTION OF LINEAR REACTION-TRANSPORT COMB EQUATION: EXTINCTION 

DYNAMICS 


The reaction term disappears from Eq. © by substitution 

P = e~ Ct P. 

The next step is the Laplace transform C[P]{t) = P(s). This yields 

sP = Dh b S{y)dlP + dpd 2 P + P 0 6(y), (7) 

and the solution can be considered as a product P(x,y,s) = h(y, s)f(x, s) The solution inside fingers is found from 
the equation in the Laplace domain 

sn(y) = d p d 2 h(y) 

and the boundary condition d y h(y)\ y -±h = 0, which yields 


h(y,s) = 


cosh 

\ h 

-1 y\)V s /d P 

cosh 

(hs/s/dp 



and n(y = 0) = 1. 

The solution is looked for in the form 


P(x,V,s) = h(y,s)f(x,s). 


that yields for the second derivative over y 

d 2 y P={-25{y)J7/dp 


sinh 

(h 

-1 y\)V s / d p 

cosh 

{hyjs/dp 



[s/d p \fi(y) | / (x, s). 


( 8 ) 
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Therefore extinct diffusion of magnesium in the balk, determined by /(a:, s), is described by equation 

/ + Po = 0 . 


Dhbdlf — 2\Jsd p tanh h\/s/d i 


(9) 


The presence of the hyperbolic tanh distinguishes two time scales/regioons. At a short time scale, when hffs/dp 1 , 
we have tanh [hy/ s/d p \ ~ 1 that yields 

Didlf - + fo = 0, (10) 

where 

Dhb Dy/aL 

* _ 2yfd~ p ~ 2Vd 

is a generalized diffusion coefficient and fo(x) = Pq{x) = Pq6(x), for simplicity. It should be admitted that the 
generalized diffusion coefficient in the x direction is independent of the finger's geometry parameters h and b. 
Performing the Fourier transform f[k) = P[f(x)], one obtains 

Po 


f(k 7 s) = 


J~s + Dik 2 

v o 


( 11 ) 


Performing the Laplace inversion of the solution m and taking into account a definition of the Mittag-LefHer function 

a 


E a ,p{z) = 


27 ri 


-dr , a, /3 > 0 , 


one obtains 




< 7+200 


e st ds 


^=Ei i (-Dik 2 ti) 
xft. 2 ’ 2 V 2 ) 


\[S + fllP yjt 


( 12 ) 


hi l* 

This solution is valid for the “short” time-scale t <C = + . Therefore, the argument of the Mittag-LefHer function 
is small, and the latter reads ism 




Dik 2 t? 

Ei i (-Dikhi) 

2 ’ 2 \ 2 / 

~ exp 

F(3/2) 


where r(3/2) = (l/2)T(l/2) = sffi/2 is the gamma function 0|. Now the Fourier inversion can be easily performed 
that yields 


h 


e- ikx e- ak2 d k = 


1" F(3/2)l 

2 

x 2 r(3/2) 

AirDit 

2 . 

exp 

D i y/t 

L 2 J 


Taking into account the reaction term, one obtains the solution for the short-time scale exstinc diffusion in the bulk 

(13) 


n/ „ Poe Ct ( yfffix 2 

P(x, y = 0,t <t 0 ) = ; exp I - 


8Dint 3 \ 2Diy/t J 

Note, that this result is independent of h and b and coincides with the subdiffusive comb limit h —> oo. 

On a long time-scale the hyperbolic tanh in Eq. d9]) has a small argument, which yields 

tanh [hyjs/d p ] ~ hyjs/d p . 

This corresponds to an equation for normal diffusion 

Dy/bjhdlf - sf + P 0 /h = 0 DyfbJh = D bh 

with a well-known solution for normal diffusion. Finally this yields the long-time scale solution in the bulk 
p(x ' 9 = M > *”>= e '“ / dx ' p ° (x,) “ p 


(14) 

(15) 


When Pq{x) = 6(x), the Green function coincides with the distribution. In Eqs. (fl3l) and (fl5l) . the transient time 
parameter j s determined from the geometry of the experiment. 
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V. REACTION FRONT PROPAGATION IN A CASE OF HIGH CONCENTRATION OF MAGNESIUM 

In a case when the concentration of magnesium is high enough, the approximation © is not valid anymore, and 
one has to take into account the second order reaction ©. Thus the reaction term reads 

C(P) = CP + CiP 2 C = p ■ k ■ (N — P 0 ) , C\ = p ■ k . (16) 

In this nonlinear case, an exact analytical treatment of Eq. © is impossible, and we apply an analytical approximation 
to find the overall velocity of the reaction-polymerization front propagation without resolving an exact shape of the 
front, namely without knowledge of the exact distribution function P(x,y,t). 

Since we are seeking for the front propagation in the x direction, the exact shape in the y direction is not important, 
and we consider the distribution as a function of the x coordinate only P\(x,t). In other words, the detail dynamics 
inside fingers is not important, and we take into account its overall contribution to diffusion in the x axis . To this 
end, the y coordinate is integrated out 




P(x, y, t)dy. 


(17) 


J —h 

Therefore, this integration in Eq. © with reaction due to Eq. (fTBl) yields an equation in the unclosed form. First of 
all, we take into account integration 


dyP(x,y,t)dy 


d v P{x,y,t) 


y=h 

= 0. 

y=—h 


This yields the equation in task 

r h 

d t Pi(x,t) = hbDdlP(x,y = 0,t) - CPi - Ci P 2 (x,y,t)dy. (18) 

J-h 

We have here two problematic terms. The first one is P{x,y = 0, f), which will be expressed by Pi(x,t). To this end 
we use the relation ©, where we pay attention on P(x , y = 0, s) = f{x , s). Integrating Eq. © over y yields 

/ h ph 

n(y,s)dy = 2P(x,y = 0,s) / h(y,s)dy. 

-h JO 


Therefore, one obtains 


P(x, y = 0, s) = y s/Adp tanh h\Js/d p 
Performing the inverse Laplace transform, one obtains from Eq. dD 

ffT+200 

ds 


Pi • 


,-sti . e ~ st 2 . e at 


(19) 


l'(J-\-lOO poo 

P(x,y = 0,t) = —— / ds dtidt 2 Pi(x,ti)TZi(t 2 ) • e~ 

2t r* Ja-ioo Jo 

/* OO poo pt 

= dti / dt 2 Pi(ti)TZi(t 2 )S(t - ti - t 2 ) = / dt'Pi(t - t')lZi(t'). 

Jo Jo Jo 


( 20 ) 


Here we take into account the causality principle, saying that P(t — ti) = 0 for t\> t and introduce the kernel TZ\ (t) 
through the Laplace inversion 


Ki{t) = C 


-l 


yjs/4dp tanh 1 (hyjs/d p J 


( 21 ) 


The second problematic term is integration of the nonlinear reaction. We present it in a form convenient for the 
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further analytical treatment. Therefore, we have the following chain of transformations: 


1 


<7+200 


N.L.R.T. = / P\x,y,t)dy= . 

J-h (ZTTl) J a _ i 

/•C7+200 


where 


e Slt dsie S2t ds 2 / P(x, y, s{)P(x, y, s 2 )dy = 


' —h 


1 


(2 tt ?:) 2 J a _ i 


e Slt dsie S2t ds 2 P(x, y = 0, Si)P(x, y = 0, s 2 ) / n{y, s 1 )h{s 2 , y)dy = 


' —h 


<><7+200 


'Jl(s 1 ,s 2 )P 1 (x,s 1 )Pi(x,s 2 )^^^- • tanh 1 h\/s 1 /d i 


(2 tt ?:) 2 

/»oo 

x e Slt e S2t ds\ ds 2 = / dTidr 2 Pi(x,Ti)Pi(x,T 2 ) 

Jo 


■ tanh 


-l 


h\ s 2 /d p 


<7+200 


( 2 tt ?:) 2 7 CT _ i 


S 2) ‘ 


y/SlS 2 

M n 


tanh 


hJsi/dp 


x tanh 


h\ s 2 /dp 


. e «i(t-n) e S2{ - t ~ T2) ds 1 ds 


P-(s 1 ,s 2 )= n(y, s 1 )h(s 2 ,y)dy. 
J-h 


Introducing new variables t — T\ = t\ and t — t 2 = t 2 yields 


dr i = — dt\ = / dt\ 


2 , 


and the same for t 2 . Therefore, we have 

nOO 

N.L.R.T. = / dt\dt 2 Pi{x, t — ti)Pi(x, t — t 2 ) 

Jo 

In what follow we will keep this form as the N.L.R.T. 
Eventually, equation for the Fi(x,t) reads 


f <7+200 


(27ri) 2 


lZ(si,s 2 )e S1 1 e S2t2 ds\ ds 2 


dtPx = hDdl f - t')dt’ - CP! - C x [ H(t\ t")Pi(t - t')Px{t - t")dt'dt" , 

J 0 J — oo 


( 22 ) 


(23) 


where two-time-point kernel TZ(t',t") is defined in the Laplace space in Eq. 


VI. HYPERBOLIC SCALING FOR THE OVERALL VELOCITY OF THE REACTION FRONT 

PROPAGATION 


To evaluate the overall velocity of the asymptotic front, we follow the hyperbolic scaling consideration, developed 
in mini and adopted in jgj for consideration of the reaction transport front propagation in comb structures. Let us 
introduce a small parameter, say e, at the derivatives with respect to time and space [l6, 17]. To this end we re-scale 

x-^x/e and t— >f/e,and Pi(x, t) —» PJ (x, t) = P\ 



. Therefore, one looks for the asymptotic solution in the form of the Green’s approximation 


Pi(x, t) = exp 


G+M) 

£ 


(24) 


The main strategy of implication of this construction is the limit e —> 0 that yields an asymptotic behavior at finite 
x and t , where we have 


exp 


G e {x,t) 

£ 


= 0 , 
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except for the condition when 


G e (x, t) = 0 . 


This equation determines the position of the reaction_spreading front, and in this limit, G(x,t) = lim^o G e (x, t) 
is accounted as the principal Hamiltonian function flti fl7| . Therefore, the Hamiltonian approach can be applied to 
calculate the propagation front velocity. In this case partial derivatives of G(x, t) with respect to time and coordinate 
have the physical senses of the Hamiltonian and the momentum: 


dG(x, t ) 
dt = “ 


d G(x, t) 
—8^ =P - 


(25) 


Now, ansatz (El) for the probability distribution function inside the bulk is inserted in Eq. (El) , where we also 
make the scaling change x —> f and t —> 

Let us start from the last term in Eq. (El) , which is the reaction term. First of all we take into account the following 
change in the upper limit of integrations 


r* ri r°° 

/ dt' =$■ lim / dt' = dt'. 

J — OO £—>0 ./ _ OC J _ go 


Then, we make the following expansion for Pf{x, | — t') = Pf (t — et '), which reads 


Pf (t — et') = exp 

Therefore, the nonlinear reaction term reads 


— G e (t-et') 

£ 


exp 


— G e (t) + t’d t G E {t) 
£ 


/ OO -I rcr+ioo 

dt!dt" ——-r / P(s', s")e~( H ~ s ^ ^ ds'ds" . 

-OO ( 27r *) Ja-ioo 


N.L.R.T. = e~ z ~ ■ / dt dt , 

(27TZ 


(26) 


where we use the first equation in (1251) . Integration over s' and s" is performed with some care, since 1Z is singular at 
s' = ±s". The main result here is that TZ is finite and N.L.R.T. ~ e -2 ^ - , that in the limit e —> 0 is of the order of 
o ) and, as will be seen in this limit, does not contribute to the finite result, since the rest of the equation is of 

the order of O ^ . 

Now we consider the the kinetic term taking the time integration with TZ\ kernel. This reads 


f £ i G £ (t — et' ,x) . G e (t) f tt ,/ 

/ 7Z \ (t ) e e dt = e ~ / lZ\{t )e~ 4 dt 

I o J o 

G e (t) ~ 

= e —TZi(H), 


where IZi(H) is defined in Eq. (12T1) . 

Finally, differentiating in the limit e — 0 where ( N.L.R.T. = 0) and taking into account that the Hamiltonian 
H and the momentum p in Eq. (1251) are independent of x and t explicitly (which leads to the absence of mixed 
derivatives), one obtains that the kinetic equation (1231) becomes a Hamilton-Jacobi equation: 


-d t G = 


h b D(d x G) 2 Pi(P) - C , 


(27) 


where G = G(x , t ) = lim e _>o G e is the action (principal Hamiltonian function) 

G(x,t) = [ \p(t)x(t) - H(p(T),x(T))]dT . 

Jo 

The rate v at which the front moves is determined at the condition G(x, t ) = 0. Together with the Hamilton equations, 
this yields 


v = x = 


dH 

dp 


_ H 
V ’ 


(28) 
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The first equation in (l28l) reflects the dispersion condition, while the second one is a result of the asymptotically free 
particle dynamics, when the action is G(x,t) = px — Ht. Taking into account x = vt, one obtains Eq. (I28fl . The 
combination of these two equations can be replaced by 


. H . H(p ) 

v = mm = mm-. 

h>o p(H) p >o p 


To proceed, we consider two time limits of the kernel/dispersion term IZi(H) in Eq. (1271) . 


(29) 


A. short time limit h^/H/d p 1 


In the short time limit, which corresponds to subdiffsion on the intermediate asymptotic times, 


tanh 


h\J H/ d p 


= 1 that yields 1Z\{H) ss sjH/4d p . Taking this into account, one obtains from Eq. (1271) 


we have 


(H + C ) 2 _ ( h b D ) 2 4 
H 4 d p P ' 

and p{H) = ■ ^ H f c . Taking logarithm from Eq. (1301) and differentiation over p , one obtains 

V hL) 4 

2 5 

H+C~H ' 


(30) 


This equation has no solution for C > 0. Therefore, for C > 0, v = min#>o = 0- This means the failure of the 
front propagation, as expected for the initial concentration of magnesium less than actin. 

For C < 0 that corresponds Pq > N, one obtains the solution H = 5|C|/3 that yields a nonzero velocity of the 
reaction front propagation at subdiffusion at the initially short times 

/5\! (2y /D 2 h 2 b \C\ 3 \i /5\! /2\ h /DaL\C\ 3 \z 

^ V3/ \3/ V 4 J p ) ~ V3/ ' V3/ V 4 d ) ' 

Again, we admit that this result is independent of the finger geometry and correspondingly, has the comb limit h —> oo. 


B. large time limit h^/H/d p <C 1 


In the opposite case of the large time asymptotics, tanh 
seen from Eqs. m and (1151) . we have from Eq. (1271) 


h\J H/ d p 


hyfH/dp, and TZi(H) 


jr. Therefore, as 


H = —^j—p 2 — C , 


and p = \J2(H + C)/Dbh■ This immediately yields 


(31) 


v = min 

H> 0 


HD 


bh 


y/2(H + C)i 


= 0 


for H = 0, as expected for C > 0. 

The situation changes dramatically for C < 0 (Pq > N ). In this case, the equation 


dH _ H 
dp p 


yields 


2Dp = Dp + J-H,. 

P 
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Here we scale the diffusion coefficient D = . This yields the solution for p = y/\C\/D. Therefore, substitution in 

the velocity equation yields 


v = nun ■ 

p>0 


H(p) Dp 2 + \C\ 


= 2JD\C\. 


p p 

This is a well known Fisher-KPP result for the reaction front velocity with the scaled diffusivity D. 


VII. CONCLUSION 

In the paper we considered anomalous transport and reaction dynamics by providing the theoretical grounds for 
the possible experimental realization of actin polymerization in comb-like geometry. Different regimes for the reaction 
of polymerization have been considered, and depending on the concentration of reagents (magnesium and actin), we 
recovered both the failure of reaction front propagation and a finite speed corresponding to the Fisher-KPP long time 
asymptotic regime. Finally, we repeat that our analytical results on reaction propagation can help to guide the design 
of microfluidic devices but also can lead to real experimental tests of anomalous diffusion and reaction dynamics. 
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